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Abstract. Most of us must have been fascinated by the eye catching displays of 
C^ ■ collectively moving animals. Schools of fish can move in a rather orderly fashion 

C*h. and then change direction amazingly abruptly. There are a huge number of further 

examples both from the living and the non-living world for phenomena during which the 
many interacting, permanently moving units seem to arrive at a common behavioural 
pattern taking place in a short time. As a paradigm of this type of phenomena we 
O ' consider the problem of how birds arrive at a decision resulting in their synchronized 

c/3 . landing. We introduce a simple model to interpret this process. Collective motion 

<*V prior to landing is modelled using a simple self-propelled particle (SPP) system with 

a new kind of boundary condition, while the tendency and the sudden propagation 

of the intention of landing is introduced through rules analogous to the random field 

Ising model in an external field. We show that our approach is capable of capturing 

^. ■ the most relevant features of collective decision making in a system of units with a 

C*~) , variance of individual intentions and being under an increasing level of pressure to 

^jj ' switch states. We find that as a function of the few parameters of our model the 

t-4- . collective switching from the flying to the landing state is indeed much sharper than 

r * the distribution of the individual landing intentions. The transition is accompanied by 

f^*) ■ a number of interesting features discussed in this report. 

o' 



PACS numbers: 87.15Zg, 87.191o 
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1. Introduction 

In recent years, there has been a growing interest in the topic of consensus decision 
making in several fields including animal behaviour [1], social sciences [2] and control 
theory [3]. In this paper we use the paradigm of animal groups but we would like to 
stress already at this point that there is a considerable analogy between the decision- 
making processes in very different systems and we expect that our simple model should 
be applicable to a wider range of systems. Our object of study is the phenomenon during 
which the animals moving in groups, seemingly make unanimous decisions on the choice 
and time of performing activities even in the absence of global leaders. It is expected 
that because of the heterogeneity in the attributes like the age, sex and social status 
of the animals in a group or the differences in their perception of external stimuli there 
will be differences in the motivations of the members at the time of making choices. 
However, in spite of the differences, synchronization is seen to take place. Examples 
include takeoff of swarm of honeybees from nest sites [I], activity synchronization in 
sheep [5], collective movement in monkeys [6], group departures of domestic geese [7] 
and departure of ants from a feeding site [8] . 

The study of collective opinion formation in physics has mainly focused on human 
societies [9]. Parallels have been drawn between phase transitions in magnets and 
the opinion dynamics in populations. Alternatively, the subject of collective motion 
of animal groups has been extensively modelled [TOj [UJ. Here point-like particles, 
representing animals, move with constant speed, each tending to align with its immediate 
neighbours, for low noise levels, giving rise to a globally ordered state replicating the 
motion of flocks where all animals move in the same direction. 

In this report we model the process of landing of bird flocks performing foraging 
flights as a typical example of collective decision making. We view this phenomena 
as a shift of the average opinion of the flock from that of continuing horizontal flight 
about some preferred altitude to that of descent towards the surface below. We model 
the fact that birds land almost synchronously even in the presence of heterogeneity 
in motivations, along the lines of the random field Ising model (RFIM) [T2l [13]. Our 
approach is similar to that used in [TJ] where the authors analyze (among other cases) 
the manner in which people in concert halls, initially joining the applause, stop clapping 
almost at the same time. In a very recent model [15] for collective landing, birds are 
assumed to move under the action of different social forces. In addition, the internal state 
of each bird is characterized by a continuous variable called landing intent such that the 
internal state of each bird is directly coupled to the internal state of its neighbours. In 
contrast we allow motivation of individual birds to be influenced by only the observable 
variables of their neighbours, like position and velocity. In addition as a part of the model 
we introduce a method which enables us to work in open boundary conditions where 
the birds in flight move cohesively and do not spread out in the horizontal plane even in 
the presence of perturbations. With this model we investigate the level of synchronicity 
in collective landing across time and space and explore the nature of fluctuations close 
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to the point in time when landing occurs. 

2. Model 

Below we describe the two aspects of the model, viz. collective decision making 
and collective motion in detail. The birds in our model are represented by particles 
characterized by the dynamical variables of position (x,y,z) and velocity (v x ,v y ,v z ) in 
three dimensional space. The only natural boundary that we allow in the system is the 
z = plane which we consider to be the landing surface. The z-coordinate of a flying 
bird gives its height above the ground. We assign to each particle i a variable ttj such 
that we refer to a bird in the flying state as active (o^ = 1) and a bird in the landed 
state as inactive (a, = 0). This variable is updated at discrete times t according to the 
condition: 

, . f 1 if zAt- At) > ... 

aAt) = < *) A ( A . (1) 

w 10 i£zi(t-At)<Oarai(t-At) = 0, 

where At is the time increment. Once a bird i has landed we keep its position unchanged 
till the end of the simulation. Also, from here on we use the words "bird" and "particle" 
synonymously. 

First we focus on the variables x^ 1 = (xj, yi) and v^ 1 = (vf, vf) governing the cruising 
motion of the birds parallel to the landing surface. This motion is assumed to be 
essentially decoupled from the motion in the vertical direction and is modelled in the 
spirit of [10J. A particle i is assumed to move, with some uncertainty, in the average 
direction of motion of all neighboring particles j whose separation |x"- — x^'| from i is 
less than a interaction radius R. From now on we refer to the set of particles in this 
neighbourhood as M^r with respect to a particle i such that, the inactive particles and 
the particle i are included as well |j. The expression for updating velocity at a time t is 
given by 



V!{t) =v-N 



M{if) ■ N( Y, V;(t - At)) + F B (x3(t - At)) 



a,-=l 



(2) 



where v is the magnitude of the velocity, Ai (£*) is the rotation matrix in two dimensions 
representing a random perturbation by an angle £| and N(u) denotes the unit vector 
u/|u|. Here £| is chosen with a uniform probability from the interval [—r]i:,T]7r], where 
rj is the amplitude of the noise. We describe the force-like term Fb below. We follow 
the standard version [17] of the original SPP model pH] for the updating of positions: 

x^)=x^-A£)+<(£)A£. (3) 

We mention here that it is a well-known problem with existing models of collective 
motion to maintain the cohesiveness of a moving flock in the presence of a finite amount 

| Neighbourhood could be defined in alternate ways as well. For example, Ballerini et al [16] found 
that for starlings, a better definition of the neighbourhood is based on rather the topological than the 
Euclidean distance 
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of noise r\ and with finite ranged (as R) interaction between the members. A flock tends 
to break up into smaller sub-flocks which eventually move in independent directions. 
A periodic boundary condition (PBC) is used to prevent the flock from spreading out 
perpendicular to the direction of motion. However, in the problem of landing we find that 
the imposition of a PBC would be unrealistic in the sense that there will be a possibility 
of a particle that has become inactive and hence ceased to move, to continually exchange 
information with the rest of the moving flock. Therefore, we introduce a new kind of 
"comoving boundary" condition where we imagine the motion of the particles on the 
xy-plane to be governed by a circular bounding region with radius R B and the centre 
lying on the centre of mass ( CoM) of all active particles (mass of each assumed to be 
unity). This boundary therefore translates with the flock as the CoM is put into motion 
due to collective motion of the particles. When a particle which is inside the region tries 
to leave it, the particle is subjected to an attraction proportional to the distance from 
the boundary. The direction of this attraction is along the resultant of the direction 
towards the location of the CoM and the direction of motion of the CoM itself. For a 
particle located at x" at an instant of time, the attraction Fb(x") is given by: 



F B (x" 
with 




Re 



C(x") if |x"-x"| > R B 

' \ (4) 

otherwise, 



C(x") = -N(x" - x") + /3N(v M ), (5) 

where D B is the strength of the attraction, x 11 and v" are the position and the velocity 
of the CoM, respectively, parallel to the xy-plane at that instant and (3 is a parameter. 
The CoM and its velocity are calculated with respect to the active particles. In the 
figure [1] we plot Fb when the motion of the CoM is in a direction parallel to the Y- 
axis. Similarly, we assume that the boundary adjusts itself to keep average separation 
between active particles constant. Thus, in general R B is a function of time, given by, 



where A(t) = ^2i(ii(t)/N is the fraction of the active particles at time t and p is the 
number of particles, initially active, per unit area of the xy-plane. 

Now we explain how the collective decision making process takes place in the flock 
and how it is related to the cohesion and eventual landing along the vertical direction. 
We characterize the decision of an individual bird as a binary choice problem [14J. To 
each bird i we assign an internal state variable Sj such that when Sj = 1 the bird 
continues to cruise above ground and when s, = — 1 the bird decides to land. For the 
cruising state we assume that a bird likes to fly about an altitude z on average without 
deviating too far above or below. However, we expect the elevation of a bird to fluctuate 
with time pertaining to factors like the noise in the medium and collision avoidance. 
Thus when s ; = 1 we consider the vertical motion of the particle to be essentially a 
random walk bounded by imaginary walls at z = z — 4p and z = z + 4p where Az is 
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Figure 1. A typical snapshot of the vector field Fb(x") on the xy-plane, defined 
through (TJ| for the value of /3 = 0.75. The direction of motion of the Co M is assumed 
to be along the Y-axis, indicated by the arrow at the centre of the circle. The length 
of a vector is directly proportional to the distance from the boundary. 



the thickness of the flock. In case the bird decides to land the vertical motion is directed 
towards the landing surface at z = 0. Formally we define, 

{e\-v if Si(t) = 1 and \6{t - At)\ < Az /2 

-v ■ sign[<5(t - At)] if Si(t) = 1 and \6(t - At) | > Az /2 (7) 
-v if Si(t) = -1, 

where e\ takes values ±1 with equal probability and 8(t — At) = z(t — At) — z , is the 
amount of deviation from the height zq at time t — At. 

To model the motivational differences in birds we assign to each bird i, an inherent 
switching time ti [IE] such that if the bird begins an isolated flight at time t = 0, it 
would decide to land at time t = t t . We assume that the value of tj for the bird i 
will in general depend on its energy reserves [19] and thus in general will be different 
for different birds. We choose tj's from a Gaussian distribution with a given standard 
deviation a Q . The value of o"o allows us to control the level of heterogeneity in a flock. 
Although we assign a priori values of t% to birds, due to the randomness in the choice of 
these switching times, they can be thought to be additionally include the effect of the 
environment such as quality of foraging patches over which the flocks fly. 

For a particle in the flock, it will have a finite number of particles in its 
neighbourhood and therefore will be influenced by their actions while taking a decision. 
Hence the internal state of a particle will be a function of two competing factors, (i) its 
inherent switching time and (ii) the internal state of the neighbours, which we assume 
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to be reflected by the nature of their motion. In an approach similar to that of the 
RFIM at zero temperature [14] we write down: 



Si (t) = sign|*i - t + JS AfiR {t - At) J (8) 

where J is the propensity of a particle to follow the decision of its neighbours and 
Sj\f t R is the local mean field through which the bird i senses the average decision in its 
neighbourhood M%,r- We expect that the inclination of a bird % to land will be enhanced 
by a resultant motion of birds in H%r towards the landing surface, and also by the 
number of birds that have already landed. We define S^ R at a time t as: 

S K Jt) = 1 + 2A M . R (t y yJ v ' K > R - 2(l - A„. R (t)), (9) 

where (v z {t))^ and A^ i R (t) are the average value of z-component of velocity of active 
particles and fraction of active particles, respectively, in the neighbourhood Afi 7 R at 
time t. At times t << (tj) — a , the birds on average have no inclination to land and 
therefore we expect them to maintain a height around z = z without any resultant 
vertical motion towards the landing surface. In this regime we expect for any particle 
h ( vZ )n- ~ an d a l so > 1 — A/V^ — 0, implying that Sj^ iR ~ 1. On the other hand a 
resultant vertical motion in negative z-direction would imply (v z )j^. — > —v. Thus in a 
regime when majority move towards the landing surface at z = and eventually reach 
inactive states, we expect Sjy. R — > — 1. We now move on to characterize the resulting 
collective landing as a function of different states of collective motion. We would like 
to emphasize at this point that although the detailed description of our model involves 
specifying nine expressions, the essential concept and the ingredients of the model are 
rather simple (SPP + RFIM) and the number of independent parameters is low (r\, R, 
a, J). 

3. Results 

We note that in our model, the shape of a coherently moving flock, far from the landing 
regime, is roughly cylindrical with a radius of cross section -Rb(O) = y/N/p and thickness 
Azq. We argue that from the point of view of our model for large flocks Azq will have in 
general weak or no dependence on the number of birds. We choose the thickness of the 
flock Az as the unit to measure distances (Az = 1.0). We choose the projected density 
p = 2.0, the interaction radius R = 2.0 and the preferred altitude of flying z = 10.0. 
In our simulations we take iV = 1024, unless mentioned otherwise, the order of which 
is typical for moderately large flocks [20] . Also the ratio Az /z would be typical in 
these flocks which have thickness of order 10 m and fly at altitudes of around 100 m. At 
the time t = in (jSJ) we expect the flock to be already in the steady state of collective 
motion. So before we switch on time as in (JSJ) we evolve a flock to its steady state of 
motion according to equations (E}, (EJ) and (JTj) where each particle i has Sj = 1. It is 
known (TOj [11] that the nature of the collective motion depends on the amplitude of the 
noise. We find that for sufficiently low values of 77, all the particles align their directions 



Collective decision making in cohesive flocks 7 

and the flock as a whole translates in some spontaneously chosen direction parallel to 
the xy-pl&ne. As the value of rj is increased, the tendency of alignment is gradually 
lost leading to weak or no translation. However, the actual nature of this transition is 
not known precisely enough and its detailed investigation is outside of the scope of the 
present paper. 

3.1. Time scales 

We make an observation that there are two independent time scales in the problem - 
the time scale of motion and the time scale over which the inherent switching times are 
distributed. The former is given by T v = Az /v which is the time required to move one 
unit of length (or half interaction radius R = 2). In presence of finite amount of noise 
rj, this time scale would also govern the mixing between particles. The width of the 
distribution of tj's (~ 2cr ) gives a measure of the time over which the energy reserves 
of the majority in the flock are exhausted. Therefore, when o"o < T v , there is hardly 
any motion during the collective decision making. When o"o 3> T v the flock actually 
translates (at low values of rj) over a long distance during the process. We call the ratio 
of these two time scales as A , where A = c /T v = cr v/Az , and we investigate the 
nature of landing at its different values. 

We set the time scale for choosing values of v and <To as the time increment At = 1.0. 
The value of magnitude of velocity v is kept around and below 0.1 so that change in the 
set of neighbours Ni,R of any bird i is gradual. We choose the value of the parameters 
for the comoving boundary described by equations (jl]) and ([5]) to be D B ~ 10~ 3 and 
(3 = 0.75. We also observe that parameter /3 controls the degree of asymmetry of the 
vector field C (J5j) outside the comoving boundary. For ft — > the system enters into 
a rotating phase [21], and for /3 > 1 the flock has a tendency to accumulate at the 
boundary. 

3.2. The landing process 

The role played by the coupling J in ([8]) is well known [TH [13] and has also been the 
matter of investigation. However, in our model we fix J to an optimal value for which 
we find the desired behaviour, as explained below. In the figure [2](a) we compare the 
process of landing across time for different values of the coupling. For a flock with 
inherent switching times {tj} chosen with a fixed cr Q we monitor the fraction of active 
particles Aj(t) with time t for different values of J. We observe that when J = the 
particles become inactive independent of each other at their own switching times as 
seen from the curve Ao(t) which falls smoothly with time. On choosing J = 10o"o the 
behaviour is quite opposite. A very sharp fall occurs in the curve AiQ ao (t) implying 
that there is indeed a collective decision making taking place when a large number of 
particles reach inactive states within a short period of time. However, comparing the 
two curves A (t) and A Wao (t) it appears that in the latter case, the collective landing 
occurs at point in time which is much larger than (by an amount of order <7o) than 
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A(t) 




12 r 0>) 







Figure 2. (a) The quantity A(i) denotes the fraction of birds not landed yet as the 
time t progresses. The curves Ao(t), A2 ao (t) and AiQ ao (t) correspond to values of 
coupling J = 0, 2ct and lOco, respectively. Also, S% ao {t) shows the fraction of birds 
in state Si = — 1 at time t for J = 2co- (b) A snapshot of a collective landing of a 
bird flock with N = 1024, p = 2.0, 77 = 0.2, i? = 2.0 and u = 0.01 with <t = 1000 and 
J = 2<7o- The arrowheads point in the direction of the velocity of the birds. 



the switching times of most of the birds. We find that a much lower value of J = 2o~q 
is relevant biologically. The curve A2a {t) sharply falls indicating collective behaviour 
and the fall occurs near at a time when about half of the birds have just crossed their 
inherent switching times. 

In the figure 12(a) we also plot S^o if) which is the fraction of birds which have a 
momentary intention to land at time t i.e., the fraction of particles with Sj = — 1. Unlike 
A2<T (t), the quantity SWoW is seen to fluctuate before becoming equal to the maximum 
value unity. This shows that before finally deciding to land, birds may intermittently 
change their decision. This is caused by the fluctuations in the vertical motion (JTj) of 
neighbours as well as the changing neighbourhood due to finite value of 77. However, 
these fluctuations become only important near the landing regime. We explore an aspect 
of these fluctuations later. In the figure E(b) we plot the velocity field of a typical flock 
when the majority has decided to land. 



3.3. Role of parameters 

Below we quantify the extent to which the presence of the coupling affects the 
collectiveness in landing as we vary the scaled value of the heterogeneity (A ). We 
define a measure of the sharpness of the decrease in the fraction of active particles 
A(t) (figure [2]). During the process of landing, we record the times at which individual 
particles become inactive. Then we calculate the standard deviation a of the distribution 
of these landing times. Finally, we divide o by o~ to get the normalized measure. 
The dimensionless quantity a/o"o averaged over different realizations of tj's, provides a 
measure of the extent to which the landing of the flock is coherent in time such that 
a = a implies the complete absence of collective decision making. In the figure [3(a) we 
plot o-/ctq versus 77 for different values of Ao- We make the following observations. For 
Aq — 1 the flock reaches consensus almost in the absence of any motion. As a result 
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Figure 3. (a) Plot of the ratio of the standard deviation (a) of actual landing times 
to the standard deviation (ctq) of inherent switching times as a function of noise T) for 
different values of Ao, i.e.., <Jo scaled by the scale of motion T v , as indicated by different 
symbols. Inset (A) shows variation of A i.e., a scaled by T v , with noise. Inset (B) 
shows the variation of (T/T ) which ignores the first and last 20% of the landings, (b) 
The variation of the spread of the landed flock over the xy-plane, relative to the spread 
at time t = 0, as a function of noise. 



there is no variation of o /o~q with the magnitude of noise. In this regime the value of 
a is governed by the spread of the flock in the steady state along the vertical direction 
i.e., Az . Therefore we expect that o ~ Az /v. In the inset (A) of the figure [3]^a) 
we plot the ratio A = a/(Azo/v) against r\. It is worth mentioning that the quantity 
A provides an alternate normalization of a with respect to the other independent time 
scale T v (— Az /v). We find that for Ao close to unity, A < 1 asserting that the landing is 
coherent and the separation between landing times is due to the width Az . For larger 
values of A the value of cr/a also becomes large signifying gradual loss of coherence. 
We also observe that in general cr/oo has a maximum for intermediate values of the 
noise. 

We find that if we choose to ignore the particles which become inactive either too 
early or too late, then the level of temporal coherence appears better. For a particular 
landing we measure the times £20 an d t 80 at which 20% and 80% of the particles, 
respectively, in the flock has just become inactive. Thus we define, T = t 80 — t 2 o as 
the width of the time window in which a majority i.e., 60% of the particles become 
inactive. We normalize T by the corresponding time window (To) for a flock with the 
coupling J = 0. In the inset (B) of the figure [3]^a) we plot (T/T ) versus r\ for different 
values of Ao- The average behaviour reveals a level of coherence such that two birds 
originally having their inherent switching times separated by T land within a period of 
around 7% of T . 

We investigate the degree of collect iveness as a function of the distance over which 
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a flock lands in figure |3](b). We calculate the following quantity: 



Ar=\l^pl-(^) 2 . (10) 

This quantifies the spread of the flock across the landing surface. We normalize the 
value of this spread (Ar) when all particles become inactive by the corresponding value 
(Ar ) calculated with the positions of the active particles at time t = 0. In the figure 
El^b) we plot (Ar/Aro) as function of noise varying values of Ao- For low and moderate 
values of Ao, the value of (Ar/Ar ) is of the order of unity. The landing occurs over a 
region whose dimensions are comparable to dimensions of the flock. 

We observe that in the regime of large Ao, the value of noise plays a major role in 
deciding the overall nature of landing. For very low values of r\ the flock translates over 
large distances. Due to this translation, the particles reaching inactive states early are 
unable to influence the decision of the active particles in the flock as the CoM moves 
away from these particles. In this case the landing occurs in few large clusters separated 
in space and time. In case of higher values of the noise there is weak or no translation of 
the CoM. However, the scale of the difference in inherent switching times i.e., <To being 
very slow as compared to the scale of motion, it provides the scope for a particle with 
a low value of ti to diffuse to the boundary. As a result, temporarily the influence of 
the neighbourhood is weakened. This allows the particle to switch its state and become 
inactive. Initially, the excursions of the particles to the comoving boundary becomes 
the reason behind independent landings. The shrinking of the boundary ([6]) effectively 
screens (depending upon the interaction radius R) the active particles from the inactive 
particles. For moderate values of r\ the path of the CoM becomes tortuous. However, 
the effect of segregation of inactive particles from the active particles still persists. The 
landing occurs in large number of smaller clusters. 

We further characterize the above behaviour by the following definition. After the 
landing we calculate the CoM of the particles positioned on the x|/-plane. We define 
the measure of spatial coherence \l/ as the fraction of particles that lie within a distance 
.Rb(O) from the CoM. Formally, 

1 N 

^"E^bW-K'-x"!), (ii) 

8=1 

where, x^ 1 is the position of the i-th inactive particle, x" is the position of the CoM and 
if (...) is the Heaviside step function. We expect from the above observations that \1/ is 
close to unity for a spatially coherent landing and zero in the opposite case. 

In the figure [4(a) we plot the positions of the birds on the xy-plsnae after all the 
birds have landed. The line shows the trajectory of the CoM and the circle is drawn 
with the radius -Rr(O) about the centre of mass calculated with the final positions of 
the landed birds. Majority of the birds lying within the circle shows that in such a 
case spatial coherence is retained. In the figure S(b) we plot \& as function of Ao for the 
value of 7] = 0.5. The different plots corresponds to different system sizes. The decrease 
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Figure 4. (a) The snapshot after the flock, described in the figure [2b), has landed. 
The line shows the trajectory of the centre of mass of the birds in flight. The arrows 
indicate the average direction of motion of the centre of mass at different points in 
time. The circle drawn is centred on the centre of mass of all landed birds and is 
drawn with the radius -Rb(O) of the comoving boundary, (b) The measure of spatial 
coherence $ as a function of Ao for the value of 77 = 0.5. Different symbols indicate 
the system sizes N = 512, 1024, 2048, and 4096. Inset (A) shows the plots with A 
rescaled to Xq/Xq(N). Inset (B) shows the variation of Xq(N) with system size N. 
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Figure 5. (a) The scaled values (A) of the standard deviation of the actual landing 
times as function of the scaled values (Ao) of the standard deviation of the assigned 
switching times for different system sizes (N = 512, 1024, 2048 and 4096) indicated 
by different symbols. Inset (A) shows the collapse of the previous plots as result of 
transformation of the horizontal axis from Ao to iV _a log(Ao/A ) with a = 0.1 and 
A = 32.0. Inset (B) shows the variation of the ratio a/ao with Ao- (b) The variation 
of A as a function of the interaction radius R for different values of Ao = 16, 128 and 
512 as indicated by different symbols, in flocks of size N — 1024. The value of 77 in 
both (a) and (b) is 0.5. 
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in \1/ with increase in Ao is consistent with all the previous observations. In the inset 
(A) of the figure @Jb) we obtain a collapse of the data after rescaling A by Xq(N). The 
variation of Aq with the system size N, as shown in the inset (B), is found to be Aq ~ iV 7 
where 7 ~ 0.6 (we note that the dimension of the flock varies as iV a5 ). This relation 
suggests that with the increase in the system size there is an increase in the value of Ao 
below which spatial coherence is expected. We investigate the dependence on system 
size of the temporal behaviour in the next paragraph. 

In the figure 0(a) we plot A versus A for four different system sizes for value of 
noise rj = 0.5. The figure reveals how the spread in actual landing times increases with 
the value of Ao- A data collapse in the inset (A) reveals a scaling relation of the form: 



A 



An n° 



A J 



(12) 



with a = 0.1 and A = 32.0. In general, we expect A to be function of the parameters 
like p and R. The equation ( ITS]) suggests that for a fixed value of Ao, with the increase in 
the system size there is increase in temporal coherence. We also plot the corresponding 
values of the ratio a /ctq in the inset (B). The plot shows how this ratio approaches the 
maximum value unity with the increase in the value of Ao- The reason for the rise in 
cr/cro for very low values of A is already stated with respect to the figure [3]^a). In all 
the above investigations we have assumed that the interaction radius controlling the 
average number of neighbours of any particle is R = 2.0. We now address the issue of 
having larger (or smaller) interaction radii. In the figure \Mh) we plot the variation of 
A with the value of R for different values of Ao with value of noise set at rj = 0.5. We 
find significant decrease in the spread A for large values of R. This suggests that global 
information can lead to a more coherent decision making process. 

3.4- Fluctuations 

Lastly, we focus on the nature of fluctuations during the landing process. As a result 
of the decision reversals in birds under the influence from the neighbourhood (as 
reflected by «S , 2 ff0 (i) in the figure EJ^a)) we find that birds originally moving towards 
the landing surface sometimes return back to the cruising altitude. We characterize 
these fluctuations by the maximum depth d max to which a bird makes an excursion 
before returning to an altitude greater than z — Az /2. In the figure [6ta) we plot the 
probability density P(d max ). We find that the distribution of d max is a power-law decay 
bounded by a minimum and maximum governed by the magnitude of velocity v and 
the maximum possible depth Zq, respectively. We also study the jumps that occur in 
position of the CoM (on the xy plane) of flying birds once the landing process starts. 
As is seen from the figure IU(a) the trajectory of the CoM becomes irregular near the 
end. The CoM actually jumps from a cluster which ceases to move to a cluster which 
continues moving. We construct the probability density P(A Co m) for the separation 
AcoM between the centres of mass at two consecutive time steps. We find that P(Acom) 
has a fat tail during the landing process. In the figure Mjo) we plot the tail of P(A Co m) 
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Figure 6. (a) The probability density function P(d max ) of the maximum depth d max 
of excursions below the flock by birds, before returning to the flock as a result of 
temporarily reverting from the decision to land. The parameter sets indicated by 
symbols are v = 0.05, ctq — 200 (circles) and v = 0.01, ctq = 1000 (squares). The 
solid line has a slope —1.85 and is a guide to the eye. (b) The probability density 
function P(Acom) of the distance travelled by the centre of mass of birds in flight, 
^CoM between consecutive time steps. The parameter sets indicated by symbols are 
same as in the case of (a) . The measured slopes of the power-law fits to the curves are 
around —2.2. In both (a) and (b) the value of Ao = aov/Az = 10, f] = 0.2 and the 
magnitude of the velocity v controls the lower cut-off of the power-law decays. 



after binning. The plots reveals power law decays with close by exponents. The power 
law has a cut-off which governed by -Rb(O). 



4. Conclusion 



In this paper we have investigated the interplay between the dynamics of collective 
motion and collective decision making. In particular we have presented a model for 
landing of bird flocks. We have identified the relevant scales in such a problem and 
we have worked in the coupling regime which is expected to be biologically relevant. 
It appears that higher values of the coupling delays the process of the global opinion 
shift. We have characterized the degree of the resulting collective behaviour by different 
measures. We found that - - as a function of the few parameters of our model - - the 
collective switching from the flying to the landing state is indeed much sharper than 
the distribution of the individual (inherent) landing intentions. Among the several 
interesting features of the model we proposed, we showed that if the flock is subject to 
larger fluctuations, the spatial spreading of the fully landed flock becomes smaller due 
to a better mixing of the information about the momentary landing intentions of the 
birds. We also find that there is a characteristic value of heterogeneity which scales with 
the size of the flock. Below this value of heterogeneity, flocks remain spatially coherent 
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in the process of landing and above this value, the coherence is lost. In addition we 
have formulated a boundary condition which is pertinent to models of flocking. We 
believe that the approach we proposed is also relevant for many other systems where 
abrupt behavioural changes are observed. Possible applications include various kinds of 
animal or even robotic groups in which an almost instantaneous global switching to a 
new desired state takes place. 
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